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ABSTRACT 

I present a 40 kilosecond Nuclear Spectroscopic Telescope Array (NuSTAR) observation of the re¬ 
cently identified low-luminosity X-ray binary and transitional millisecond pulsar (tMSP) candidate 
IRXS J154439.4—112820, which is associated with the high-energy y-ray source 3FGL J1544.6-1125. 
The system is detected up to ^30 keV with an extension of the same power-law spectrum and rapid 
large-amplitude variability between two flux levels observed in soft X-rays. These findings provide 
further evidence that IRXS J154439.4—112820 belongs to the same class of objects as the nearby 
bona fide tMSPs PSR J1023-1-0038 and XSS J12270-4859 and therefore almost certainly hosts a mil¬ 
lisecond pulsar accreting at low luminosities. I also examine the long-term accretion history of IRXS 
J154439.4—112820 based on archival optical, ultraviolet, X-ray, and y-ray light curves covering the 
past ^decade. Throughout this period, the source has maintained similar flux levels at all wave¬ 
lengths, which is an indication that it has not experienced prolonged episodes of a non-accreting radio 
pulsar state but may spontaneously undergo such a state transformation in the future. 

Subject headings: pulsars: general — stars: neutron — X-rays: binaries 


1. INTRODUCTION 

The recent discovery of three millisecond pulsar (MSP) 
binary systems that alternate between clearly distin¬ 
guishable rotation- and accretion-powered states have 
revealed a ne w aspect of compac t binaries containing 
neutron stars d^aoitto et al.ll20I3l iPatrnno et al.l [2?iT4 
iBassa et~ .1201411 . In their accreting states, these objects 
show rapid optical /UV variability and a bimodal X-ray 
variability pattern (iLinares et al.ll20l4[de Martino et al.l 
I20I3I iBogdanov et al.ll20I5ll —the X-ray flux switches 
rapidly between two distinct high and low “modes” that 
differ by a factor of 6—10. In addition, these systems con¬ 
stitute the only y-ray-emitting variety of low-mass X-ray 
binary; both confirmed nearby transitional MSPs, PSR 
J1023-b0038 and XSS J12270-4859, are bright Fermi 
LAT sources even when accreting. The unique com¬ 
bination of observational characteristics offers a fairly 
straightforward way to identify additional objects that 
belo ng to this class. _ 

In IBogdanov fc Haloernl (|20I5ll . we examined the X- 
ray and optical/UV properties of the bright ROSAT 
source IRXS J154439.4—112820, which is position- 
ally coincident wit h the high-energy y-ray source 
3FGL J1544.6-1125 (jStephen et al.l [201(11 iMasetti et al.l 
1201311 . The variability and spectral properties of 
this system were found to be essentially identical to 
those observed from the nearby tran sitional millisec¬ 
ond p ulsar binaries PSR J102 3-I-0038 (IBogdanov et al.l 
and XSS J12270-4859 (Ide Martino et all 1201311 . 
The available observational evidence points to IRXS 
J154439.4—112820 most probably being a transitional 
MSPs in a low-luminosity accreting state, making it only 
the fourth such system to be identified. 

Two of the transitional MSPs, XSS J12270-4859 


and PSR J1824-2452I, are currently in a dormant, 
accretion-disk-free state (IBassa et al. 1201 4 iRov et al.l 


120151: IBogdanov et al.ll2014al : iPanitto et al.ll2013|l . leav- 
ing only PSR J1023-b0038 and IRXS J154439.4-112820 
available for further detailed studies of the accreting 
state. At present, the underlying physical processes re¬ 
sponsible for the observed phenomenology, especially the 
rapid X-ray mode switching and unexpectedly bright 
high-energy y-ray emission associated with the accret¬ 
ing state, are poorly understood. In light of the transi¬ 
tory nature of the accreting state, it is therefore crucially 
important to investigate each system extensively at all 
wavelengths before a transformation to a possibly pro¬ 
longed non-accreting episode (lasting years to decades). 

On-going multiwavelength stud ies of PSR JI023 - t-0038 
in its present accreting state ([Archibald et al.l 120151 : 
iDeller et al.l 120151 : IBogdanov et al.l I2015D are revealing 
new information regarding the accretion and jet pro¬ 
duction processes that operate in transitional MSPs. In 
this paper, I present complementary results on IRXS 
J154439.4—112820 including a target of opportunity 
hard X-ray observation using the Nuclear Spectroscopic 
Telescope Array (NuSTAR), which further confirm the 
similarities of this system with transitional MSPs. I also 
investigate the broad band spectral energy distribution 
and long-term variability of the system, which further 
strenghten the classification of this system as a transi¬ 
tional MSP. This work is organized as follows. In §2, 
I summarize the NuSTAR observation and data analy¬ 
sis. In §3, I describe the X-ray spectroscopic analysis, 
while in §4 I focus on the hard X-ray variability. In §5 
I discuss the broad band spectral energy distribution of 
IRXS J154439.4—112820. I offer a discussion in §6 and 
conclusions in S7. 


^ For the third transitional MSP, PSR J1824—24521, which is 
situated in the globular cluster M28, its 7 -ray emission cannot be 
reilabily disentangled from the 11 other radio MSPs in the cluster. 


2. OBSERVATION AND DATA ANALYSIS 

The IRXS J154439.4—1128 20 sy stem was observed 
with NuSTAR (|Harrison et al.l[2013f) on 23 March 2015 
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Fig. 1.— NuSTAR FPMA+FPMB (black) and XMM-Newton 
MOSH-MOS2 (red) spectra of IRXS J154439.4—112820jointly fit¬ 
ted with an absorbed power-law model. The bottom panel shows 
the best fit residuals expressed in terms of a with error bars of size 
unity. See text for best fit parameters. 

(ObsID 90001007) in a 40 ks effective on-source exposure 
obtained through the Director’s Discretionary Program. 
NuSTAR is the first focusing hard X-ray telescope, with 
sensitivity over the 3-79 keV energy range. It consists 
of two co-aligned telescopes, which use grazing-incidence 
mirror assemblies to focus hard X-rays onto two Focal 
Plane Modules, FPMA and FPMB. 

The data were processed using the nupipeline script 
in NuSTARDAS and the images, spectra, and light 
curves were generated using the nuproducts task. IRXS 
J154439.4—112820 is the only hard X-ray source within 
the telescope field of view and is detected with high 
significance up to ^30 keV. For both the spectroscopic 
and time variability analyses, the source counts were ex¬ 
tracted from a circular region of radius 60" (which en¬ 
circles ~80% of the total s ource energy) centered on th e 
optical position reported in iBogdanov &: Haloernl (I2015D . 
The background was taken from a 60" radius source-free 
region. 

3. X-RAY SPECTROSCOPY 

I co nducted a spe ctral analysis in XSPEC version 
12.8.2 (IArnaudlll996f) . To maximize the photon statis¬ 
tics, the FPMA and FPMB spectra and response func¬ 
tions were combined using the addascaspec command in 
FTOOLS. Based on the results from the XMM-Newton 
spectroscopy presented in (|Bogdanov fc HalDernll2015ll . 
the NuSTAR FPMA and FPMB spectra were fitted 
with an absorbed power-law assuming the tbabs inter- 
stellar absorption m odel and chemical abundances from 
I Wilms et al.l (120001) . This produces a best fit with hydro¬ 
gen column density of A^h < 1.3 x 10^^ cm“^, spectral 
photon index F = 1.68 ± 0.06, an unabsorbed 3 — 79 
keV flux of Fx = (6.6 ± 0.4) x 10“^^ erg cm“^ s“^, and 
= 1.14 for 156 degrees of freedom. Owing to the 
hard response of NuSTAR, the fit is insensitive to the 
value of N and only a crude upper limit is obtained. 
To further refine the spectral parameters, I also con¬ 
ducted a joint fit to the XMM-Newton and NuSTAR 
data. The resulting best fit with an absorbed power-law 
gives A^h = (1.7±0.1) x lO^^ cm-^, F = 1.67±0.02, and 


Fx = (8.4±0.2) X 10“^^ erg cm“^ s“^ (0.3-79 keV), with 
xt = 0.96 for 471 degrees of freedom. All uncertainties 
quoted are at a 90% confidence level. No adjustment to 
the cross-normalizations between the two data sets was 
required. As apparent from Figure 1, the spectral con¬ 
tinuum shows no evidence for a break or turn-over within 
the NuSTAR band. 

I also explored the possibility of spectral variability 
by computing the hardness ratio variations as a func¬ 
tion of count rate. There are no statistically signifi¬ 
cant differences in hardness between the high and low 
flux modes, which is consistent with the result from the 
NuSTAR data of PSR J1023±0038 in its accreting state 
(|Tendulkar et al.ll2014ll . 

4. HARD X-RAY VARIABILITY 

A total background-subtracted hard X-ray light curve 
was generated by combini ng the emission from FPMA 
and FPMB. As shown in IBogdanov fc Halpernl (I2015f) . 
the soft X-ray emission exhibits rapid switching between 
two flux modes, with ingress and egress durations lasting 
of order 10 s. Despite the much lower count rate com¬ 
pared to the 0.3-10 keV XMM-Newton light curve, the 
rapid variability is still clearly apparent. The 3-79 keV 
light curve of IRXS J154439.4—112820, binned at 100 s 
resolution, is shown in Figure 2. While the source spends 
most of the time at a high flux level (~0.1 counts s“^), 
multiple instances of low flux levels (reaching down to 
'^O.Ol counts s“^) are evident. At the 100 s bin resolu¬ 
tion some short-lived low modes are likely averaged with 
the high modes occuring within the same time bin. A 
Kuiper’s test for time variability on the unbinned time 
series, corrected for good time intervals by removing gaps 
in exposure, gives a 9 x 10“^® (w7.9(t) probability that a 
distribution with a constant count rate would exhibit the 
observed level of non-uniformity. Due to the relatively 
low source count rate, the unknown orbital period, as 
well as known calibration issues with the NuSTAR on¬ 
board clock, I have not attempted a periodicity search. 

5. LONG-TERM BROADBAND BEHAVIOR 

PSR J1824-2452I in M28 was observed to be in a very 
faint {Lx ~ 10^^ erg s“^) non-accreting state in 2002 
(|Becker et al.ll200^ and a low-lu minosity (Lx ^ 3 x 10^ ^ 
erg s“^) accreting state in 2008 l|Bogdanov et al.l [201111 . 
The optical bright ness between these tw o states differed 
by ^2 magnitudes (IPallanca et al.ll2013ll . As reported in 
iPatruno et al.l (|20I4fl . the re-apoearance of the accretion 
disk in PSR JI023±0038 was accompanied by a ~3.5 
magnitude increase in UV brightness and a ^20-fold rise 
in X-ray luminosity. Similarly, in the state transition 
of XSS J12270-4859 (iBassa et al.ll2014ll . the transforma¬ 
tion from an accreting state to a radio pulsar state was 
signaled by a decrease in X-ray luminosity by a factor of 
'^lO and decline in UV and optical brightness by ^^3 and 
~2 magnitudes, respectively. For PSR JI023±0038, the 
return of the accretion disk also coincided with a factor of 
^5 enhancement in high-energy y-ray luminosity as seen 
by Fermi LAT (|StaPDers et all [201^ . The disappear¬ 
ance of the accretion disk in XSS JI2270-48 59 was asso- 
ciated with a ^2-fold decline in y-ray flux (lYi fc Wand 

As IRXS J154439.4—112820 closely resembles these 
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Fig. 2.— NuSTAR FPMA+FPMB 3-79 keV background-subtracted light curve of IRXS J154439.4—112820binned at a 100 s resolution. 
The regular ~3 ks gaps are due to Earth occultation. Note the occasional short-lived rapid drops in count rate, especially towards the end 
of the observation. 


transitional MSP systems, it may have also experienced 
transitions to and from a disk-free state in the past. 
Examination of archival data is particularly valuable in 
this regar d, considering t hat both PSR J1824-2452I (in 
2008-2009 iPallanca et al.|l2013l: iL inares et al.ll20l4l and 
XSS J12270-4859 (in 2012 lBassa et al.ll2014ll were uncov¬ 
ered retroactively to have undergone state transitions in 
archived data sets. Prompted by this possibility, I have 
investigated the optical/UV, X-ray, and y-ray variabil¬ 
ity of IRXS J154439.4—112820 on time scales of years. 
For this purpose, I use publicly avail able data from th e 
Catalina Real-time Transient Survey ([Drake et al.ll2009ll . 
Swift XRT and UVOT, and Fermi LAT data. Figure 3 
shows the lightcurves from these data sets, which are 
described below. 

5.1. Catalina Real-time Transient Survey 

The Catalina Real-time Transient Survey (CRTS 
iDrake et al.1120091 1 acquires unfiltered exposures of the 
sky away from the Galactic plane that usually reach 
equivalent V filter of ^19-20. The CRTS archive 
includes 324 observations of IRXS J154439.4—112820 
spanning from late 2005 to late 2013. Removing points 
with magnitude errors greater than 0.2 results in 260 in¬ 


dividual exposures. 

5.2. Swift UVOT and XRT 

The error ellipse of 3FGL J1544.6-1125 has been tar¬ 
geted with Swift once in 2006 and 2007, twice in 2012, 
and on six occasions in 2015. 

The Swift UVOT data include four exposures with 
the UVWl, two with the UVWl, three with the UM2, 
and three with the U filter. The photometric measure¬ 
ments and associated uncertainties were extracted with 
the uvotsource task in FTOOLS, using a source extrac¬ 
tion reg ion of radius 5" centered on the optical position 
given in iBogdanov fc Haloernl (|2015| 1 . 

The 12 Swift XRT observations were all obtained in 
photon counting (PG) mode. The 0.3-10 keV light 
curves were extracted using the xrtpipeline command 
in FTOOLS. The counts were grouped adaptively to en¬ 
sure that each time bin contained at least 20 counts. 

5.3. Fermi LAT 

To investigate the high-energy y-ray variability of 
IRXS J154439.4—112820 on long timescales, I have re¬ 
trieved all Pass 8 Fermi Large Area Telescope (LAT) 
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Fig. 3.— Long-term light curves of IRXS J154439.4—112820 from 
the Catalina Real-time Transient Survey (CRTS), Swift UVOT, 
Swift XRT 0.3-10 keV, and Fermi LAT 0.1-100 GeV (from top to 
bottom, respectively. The different symbols for the Swift UVOT 
data correspond to the various filters: U (solid circles), UVWl 
(open squares), UVM2 (open circlesO, and UVW2 (solid triangles). 

data within 20 ° of the source from the start of the mis- 
sion on 2008 August 4 through 2015 June 7. The 7 - 
ray light curve was obtained by dividing the data set 
in 60 day intervals and carrying out a binned likelihood 
analysis for each time segment using the Fermi Science 
Toolf0 vlOrOpS. Following the analysis procedures rec¬ 
ommended by the Fermi Science Support Center, I first 
generated counts, exposure, and source maps, livetime 
cube, and an input source model based o n the source 
inform ation provided in the 3FGL catalog (lAcero et alJ 
I2OI5II . In the likelihood fitting with the gtlike task, 
the contribution from all sources within 20° of IRXS 
J154439.4—112820, as well as from the Galactic and ex- 
tragalactic diffuse emission, was taken into account. The 
best fit values and uncertainties of the flux of 3FGL 
J1544.6-1125 were used to construct a light curve in the 
0.1-100 GeV band. 

Although the optical, UV and X-ray data have gaps in 
coverage spanning several months to years, when com¬ 
bined with the 7 -ray data, they provide useful constraints 
on the long-term behavior of IRXS J154439.4—112820 
within the past decade. For instance, it is apparent that 
there have been no intervals spanning months to years in 
which the system has been in a low-flux, non-accreting 
state. The Fermi LAT 0.1-100 GeV light curve shows 
no appreciable deviations from the mean photon flux of 
1.8 X 10“^ cm“^ s“^. This is consistent with variabil¬ 
ity index of 47.5 reported for 3FGL J1544.6-1125 in the 
Fermi LAT 4-year catalog (lAcero et al.ll201^ . which in¬ 
dicates no statistically significant flux changes over the 
course of the mission. 

On the other hand, the GRTS and the most re- 


^ Available at http://fermi.gsfc .nasa.gov/ssc/data/analysis 


cent Swift UVOT UVWl data show instances of en¬ 
hancement in brightness by ^1 — 1.5 magnitudes com¬ 
pared to the typical value. These events may cor¬ 
respond to o ccasional flaring episo des seen from PSR 
.11023 +6638 (IRogdanov et al I [ 2 m^ and XSS J12270- 
4859 (|de Martino et al.l 1201, 3fl in the optical, UV, and 
X-rays. The Swift XRT 0.3-10 keV exposures exhibit 
large-amplitude variations that span a factor of ~10 in 
count rate. The corresponding range of fluxes are in good 
agreement with the high and low mode levels o bserved 
with XMM-Newton (jBogdanov fc Halpernlf^oT^ . 

6 . DISGUSSION 

The hard X-ray properties of the IRXS 
J154439.4—112820 system as revealed by NuSTAR 
are virtually identical to thos e of the confirmed tran- 
sitional MSPs. For example, Ide Martino et all (j201(l 
[ 2 ?fn investigated the broad-band X-ray properties 
of XSS J12270-4859 during its accreting state using 
data from XMM-Newton, RXTE, and INTEGRAL. 
The source spectrum in the 0.1-100 keV can be fitted 
by a single power-law with index F ss 1.7. In the 
RXTE data, XSS J12270-4859 showed the same mode 
switching and f l aring behavior found in soft X-rays. 

I Tendulkar et al.l (|2014ll presented NuSTAR observa¬ 
tions of PSR J1023+0038 in both the radio pulsar 
and accreting X-ray binary states. The average hard 
X-ray emission in the more luminous accreting state is 
well-described a power-law with photon index F = 1.66 
with a luminosity of 6 x 10^^ erg s“^ (3-79 keV). Rapid 
variability between two clearly separated flux levels was 
also found. 

The available optical/UV, soft and hard X-ray and 
high-energy 7 -ray data for IRXS J154439.4—112820 al¬ 
lows an examination of the broad-band emission spec¬ 
trum of the system. Figure 4 shows the optical to GeV 
7 -ray range for IRXS J154439.4—112820. For an illustra¬ 
tive comparison, the multi-wavelength data set of PSR 
J1023+0038 in its accreting state is also shown. The 
overall similarities in the shape of the broadband spectral 
energy distributions offers an additional indicator that 
IRXS J154439.4—112820 belongs to the same class of ob¬ 
jects as PSR J1023+0038, XSS J12270-4859, and PSR 
J1824-2452L In light of this, the same interpretations 
used for these objects can be invoked to explain the dif¬ 
ferent emission components of IRXS J154439.4—112820. 
The optical flux is likely a combination of emission from 
the companion and the accretion disk, with the latter 
dominating the emission in the UV. The X-ray emission 
presumably originates in the inner disk, close to or at 
the surface of the neutron star. The detection of coher¬ 
ent X-ray pulsations dur ing the high mode X- ray emis¬ 
sion in PSR J1023-I-0038 (jArchibald et al.l[2015l l and XSS 
J12270-4859 (Pa pitto et al. 2015) ari d the absence of 
radio pulsations (iStapoers et al.l 12014 iBogdanov et al.l 
I2OI5II suggests that at least during these intervals active 
accretion onto the stellar surface is taking place and the 
radio pulsar mechanism is quenched. The cause of the 
aperiodic rapid drops to a low flux mode is less clear as 
it ma -v correspond to changes to a propeller ejection pro¬ 
cess (flllarionov fc SunvaevI 197^ or a dead/trapped disk 
regime (iSiuniaev fcShakural 119771 : iSpruit fc TaamI 119931 : 
ID Angelo fc Spruitll2010l . 1201211 . Moreover, it is not clear 
/soft^g^j^gj. ^;];^g inflow of matter is via a “classical” thin 
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Fig. 4.— The broadband spectrum of IRXS J154439.4—112820(black data points) covering the optical to 100 GeV 7 -rays range, showing 
the Fermi LAT (open stars), NuSTAR (open squares), XMM-Newton EPIC (solid circles), Swift UVOT (solid squares) an d XMM-Newton 
OM B filter (solid triangles). The light grey points correspond to the spectral energy distribution of PSR 11023-^0038 (see lBogdanov et al.l 
120151. for details). 



disk o r a radiatively inefficient accretion flow (I Rees et alJ 

[HH). 

The persistent 7 -rays from this system possibly arise 
due up-scattering of X-rays at the boundary between 
the accretion flow and the neutron star magnetosphere. 
Two general varieties of models have been put forth to 
explain the observed 7 -rays, which differ principally in 
the assumption of an active or inactive rotation-powered 
pulsar wind. In the active pulsar wind scenario, the X- 
ray and 7 -ray emission arise due to an intrabinary shock 
produced at the interface of th e accretion disk and th e 
pulsar wind (jTakata et al.l[2?)l^ iCoti Zelati et al.ll2ni4ll . 
The X-rays are generated via synchrotron radiation in 
the shock, while the 7 -rays are the product of inverse 
Compton scattering of the photons emanating from the 
acc retion disk off the Pulsar wind particles. As discussed 
by (jPaoitto fc Torreal2015h . for the case of an inactive 
(“quenched”) pulsar wind, the system is assumed to be 
in the so-called propeller regime in which most of the in¬ 
flowing material is expelled by the rapidly rotating pul¬ 
sar magnetosphere. The X-ray luminosity is then due 
to synchrotron plus disk emission, with the flux mod- 
ing arising from rapid switching between accretion onto 
the star (the high mode) and propeller ejection (the low 
mode). The 7 -ray emission is produced by synchrotron 
self-Compton emission buy accelerated particles in the 
turbulent boundary between the pulsar magnetosphere 
and the accretion disk. 

While IRXS J154439.4—112820 appears to be a close 
analog of the confirmed transitional MSPs, some notable 
differences are also apparent. In particular, the 7 -ray 
flux of IRXS J154439.4—112820 is ~3 higher that that 


of PSR J1023-I-0038, whereas the mean X-ray flux is ~3 
times lower. While the lower X-ray flux can be attributed 
to a greater distance (assuming a comparable intrinsic 
luminosity), the relative differences in X-ray to 7 -ray 
flux may arise either due to higher efficiency of the 7 - 
ray production mechanism or due to a more favorable 
viewing angle if the radiation pattern is not isotropic. In 
addition, for PSR J1023-I-0038, the high and low mode 
X-ray fluxes differ by a factor of Ri 6 , while for IRXS 
J154439.4—112820 the flux varies by a factor of RilO. 
Assuming that during the low X-ray mode, the accre¬ 
tion flow is cleared out of the pulsar magnetosphere, this 
difference could be ascribed to a larger accretion disk 
truncation radius in IRXS J154439.4—112820. If this 
radius coincides with the pulsar light cylinder (cP/27r) 
it implies a more slowly spinning neutron star, although 
other factors like a less massive neutron star or a dif¬ 
ferent magnetic field configuration may play a role in 
determining the high/low flux ratio. 

7. CONCLUSIONS 

I have presented the first hard X-ray observation of 
the low-luminosity X-ray binary and Fermi LAT source, 
IRXS J154439.4-112820 / 3GL J1544.6-1125. The bi¬ 
nary is detected at high significance up to ~30 keV with 
a continuation of the power-law spectrum observed in 
soft X-rays. The hard X-rays also exhibit variability on 
timescales of 10-100 seconds by a factor of RilO in flux 
as seen in the 0.3-10 keV range with XMM-Newton. The 
spectral energy distribution of IRXS J154439.4—112820 
covering the optical through GeV range is qualitatively 
very similar to that of PSR J1023-I-0038 (as well as 
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XSS J12270-4859; see Ide Martino eJ^ [Ml . 

These findings further strengthen the argument that 
this system belongs to the same class of objects as 
PSR J1023+0038, XSS J12270-4859 and PSR J1824- 
241, namely a compact binary hosting a MSP that is 
presently accreting at low luminosities. 

I find no clear evidence in archival data for prolonged 
low flux intervals at any wavelength implying that within 
the past ~decade IRXS J154439.4—112820 has most 
likely not experienced a transformation to a disk-free 
state. Nevertheless, as the strongest transitional MSP 
candidate, IRXS J154439.4—112820 warrants dedicated 
monitoring to identify behavior that may signal the dis¬ 
appearance of its accretion disk and the activation of the 


rotation-powered millisecond pulsar. 
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